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We have evaluated the energies required to twist carbon 
nanotubes (NTs), and investigated the eflects of these distor- 
tions on their electronic structure and electrical properties. 
The computed distortion energies are high, indicating that it 
is unlikely that extensive twisting is the result of thermal ex- 
citation. Twisting strongly affects the electronic structure of 
NTs. Normally metallic armchair (n, n) NTs develop a band- 
gap which initially scales linearly with twisting angle and then 
reaches a constant value. This saturation is associated with a 
structural transition to a flattened helical structure. The val- 
ues of the twisting energy and of the band-gap are strongly 
affected by allowing structural relaxation in the twisted struc- 
tures. Finally, we have used the Landauer-Biittiker formalism 
to calculate the electrical transport of the metal-NT-metal 
system as a function of the NT distortion. 

Research on the electronic properties of carbon nan- 
otubes (NTs) has been an especially active field in recent 
years. The reasons for this are twofold: on one hand they 
are viewed as potential building blocks of nano-sized elec- 
tronic devices On the other hand, they are ideal 
materials to study electrical transport phenomena in low- 
dimensional systems, and to test the theoretical models 
proposed to explain such phenomena. NTs are usually 
thought of as a graphene sheet rolled up to form a com- 
pact cylinder with a radius that can be as small as a 
few Angstroms It is obvious, however, that this pic- 
ture of a NT as a straight, geometrically and atomically 
perfect tube is somewhat oversimplified. Images of NTs 
quite often reveal structural deformations such as bent 
twisted ||], or collapsed tubes. These deforma- 
tions may develop during growth, deposition and pro- 
cessing, or following an interaction with surface features 
such as electrodes, or other NTs. An investigation of the 
electronic structure of weakly distorted nanotubes has 
been performed using a low energy field theory descrip- 
tion !§. Here we aim to study the distortion energies, 
electronic properties, and electrical transport of carbon 
nanotubes under a wide range of twisting distortions. 

The electronic properties of carbon nanotubes are usu- 
ally discussed in terms of the tt electronic structure of the 
2D-graphene sheet. The mapping of the graphene sheet 
onto a cylindrical surface is specified by a superlattice 
translational vector T = nTi + mT2 where Ti = a(l, 0) 
and T2 = a(l/2, •\/3/2), and a is the length of the primi- 
tive translation vector of the graphite lattice. The pair of 
indices (n, m) describes how the NT is wrapped to form 
the cylinder, and determines its electrical properties. An 
(n, n), so called "armchair" tube, has only two linearly 



independent Fermi points iif, where K = (kpjO). It 
is a metal with two bands intersecting the Fermi energy 
|l0|-p^. Thus, when an {n,n) NT is atomically perfect 
and undistorted, it is expected to exhibit at low bias a 
quantized conductance equal to 4e^//i. However, even in 
the case of a perfect NT, a reduction of the conductance 
may be observed due to the non perfect transmission at 
the NT-metal electrode contacts. 

The low energy theory ||,^ used previously to describe 
the electronic structure of distorted carbon NTs considers 
only the 7r-orbitals, and can be appropriately described 
by a tight-binding model in which deformations are taken 
into account simply by modifying the hopping ampli- 
tudes at the deformation sites. This analysis shows that a 
smooth bending of the tube along its axis of rotation has 
little influence on its electronic properties. On the other 
hand, it predicts that a twisting deformation can lead 
to substantial modification of the electronic structure. 
The low energy theory is applicable to small deforma- 
tions. Large deformations induce the mixing of orbitals, 
so that describing the electronic structure changes simply 
by modifying the hopping amplitudes is no longer ade- 
quate. In previous work |13| , we computed the electronic 
properties of bent NTs. We found that there is a max- 
imal bending angle above which the low energy theory 
breaks down and the transport properties are substan- 
tially changed. Here we extend this work by studying 
the effect of twisting distortions over a wide range of 
twisting angles using a model that takes into account of 
both the s and p carbon orbitals, as well as allowing the 
geometrical relaxation of the distorted structures. 

The nanotube model used in the computations con- 
tains 948 carbon atoms arranged in an armchair (6, 6) 
structure. The energetics of the deformations were de- 
termined with molecular mechanics using the TINKER 
program |Q with a modified MM3 force- field |15|. Cal- 
culations were also performed using a tight-binding den- 
sity functional theory (TB-DFT) method developed by 
Porezag et al. ||l^ . In both types of calculations the dan- 
gling bonds at the ends of the tubes were saturated by 
hydrogen atoms. Distorted structures were generated by 
two different schemes. In the first, continuous twisting, 
the circular plane sections of carbon atoms along the tube 
axis were rotated sequentially by an (additive) specified 
angle. In the second, alternate twisting, a (non-additive) 
rotation was performed on every second section. The 
molecular mechanics calculations showed that it is less 
energetically demanding to continuously twist the tube 
than to twist it alternatively. Therefore, in the rest of the 
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paper, we will focus our attention on continuous defor- 
mations. The electronic structures of twisted NTs were 
determined using the extended Hiickel method (EHM) 
[ prt - EHM gives results similar to those obtained on ex- 
tended NTs with more sophisticated methods The 
energy band-gap was obtained as the difference between 
the HOMO and LUMO energies. 

To evaluate the transport properties, the two ends 
(with dangling bonds) of the NT were bonded to gold 
electrodes. The computation of the transmission function 
was then carried out by using the Landauer-Biittiker for- 
malism as described in detail in ref. In this scheme, 
the retarded Green's function for the system is given by: 



where H^t and Snt are the Hamiltonian and the over- 
lap matrices for the carbon nanotube, and Ei 2 are self- 
energies describing the interaction of the NT with the 
metal contacts. These are given by: 

^3 = T]9j-rj (2) 

with Tj describing the coupling between the NT and con- 
tact j and gj the Green's function of the metal used for 
the contacts. The total transmission function from one 
contact to the other is then given by: 

T{E) = TrpaGTiG^] (3) 

with Tj = — sj). 

In Figure 1, we plot the potential energy for a contin- 
uous twisting deformation against the twist angle for a 
constrained (6, 6) structure, and for a relaxed one where 
only the first two carbon sections at each end of the tube 
were held fixed. In the inset, we compare the energies 
obtained with molecular mechanics for small angle twist- 
ing of the relaxed structures to the energies determined 
with the TB-DFT method. A number of observations 
can be made. First, the energy required for even a mod- 
erate twisting of the tube is very large. At low twisting 
angles (inset) a quadratic law is obeyed. The twisting 
energies calculated using the MM3 potentials are lower 
than those computed using TB-DFT (inset). Therefore, 
in the following we will consider the MM3 energies as 
lower limits to the energetics of twisting deformations. 
Figure 1 shows clearly that structural relaxation has a 
very strong effect on the potential energy of deformation, 
especially for large twisting angles where the relaxation 
energy is computed to be of about the same order of mag- 
nitude as the final energy. Furthermore, the energy pro- 
file of the relaxed structure deviates from the quadratic 
law (dot-dashed line) for 87- > 14°nm^^. This deviation 
corresponds to the onset of a structural transformation 
where the tube flattens and takes on a helical shape (such 
as in Figure 2 [c,d]). A similar transition has been ob- 
served in molecular dynamics simulations [pO[. In order 



to determine the effect of the diameter on the collapsing 
of the nanotubular structures, we have also computed 
MM3 energies of twisted (10, 10) models. Taking into 
account the differences in the number of carbon atoms 
involved ((10, 10)/(6, 6) = 1.67), and the NT diameters 
('/'(io,io)/'/'(6,6) — 1-67), we find that the energy needed 
to twist a (10, 10) NT is similar to that needed to twist a 
(6,6) NT. The collapse of the (10, 10) tubular structure 
into the helical shape occurs at a lower twisting angle 
than the (6, 6) structure, i.e. when 8t > 10°nm~"'^. How- 
ever, due to its larger diameter, the relative displacement 
of the carbon atoms during the tube to helix transition 
for a (10, 10) NT is larger than in the case of a (6, 6) NT. 
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FIG. 1. Molecular mechanics calculation of the energy re- 
quired to twist a (6,6) armchair NT as a function of twist 
angle, for a constrained (open circles) and a relaxed geome- 
try (closed circles). The dot-dashed line corresponds to the 
quadratic law determined from the low twisting angle region 
of the relaxed NT. The inset shows a comparison of the MM3 
and TB-DFT energies for small twisting angle. 

Our calculations also showed that the twisting ener- 
gies scale nearly linearly with the length of the nanotube. 
Therefore, we can extrapolate and obtain the twisting en- 
ergies of longer nanotube segments, and use these ener- 
gies to evaluate their Boltzmann distributions as a func- 
tion of temperature. The relative population of twisted 
(6, 6) structures can be described by: 

^^(eT)=exp[-a.^] (4) 

where, a is a constant equal to 580 nm-degrec^^-K, L 
is the nanotube length in nm, Ot is the twisting an- 
gle (degree-nm~^), and T is the temperature. From this 
equation, it is clear that large twists cannot be generated 
by thermal excitation. For example, to obtain by ther- 
mal excitation a relatively small population (1%) of a 1 
/im long (6,6) nanotube twisted by 9.8°nm~^, which is 
equivalent to the twist deformation obtained by Clauss 
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et al. from STM images of nanotube ropes 1^, it would 
require a temperature larger than 10^ K! We conclude 
then that large scale nanotube twists must be the re- 
sult of either mechanical interactions taking place during 
the deposition or manipulation of the nanotube sample, 
or they are introduced during the high energy growth 
processes and frozen in place by shear forces due to the 
interaction with the substrate and other tubes. Local 
twisting associated with tube collapse is, however, less 
energetically demanding and is encountered quite often 
in AFM images. 




abed 

FIG. 2. Relaxed NT structures obtained with force-field 
MM3 energy minimization on (a), 1 (b), 2 (c) and 
3 "section"^ (d) twisted models. 

We now consider the electronic structure changes in- 
duced by twisting the NTs. Figure 3 shows the variation 
of the band-gap value of a twisted (6, 6) nanotube as a 
function of the twisting angle, for both constrained and 
relaxed NT structures. We find that for the constrained 
structure, the band-gap increases linearly with twisting 
angle up to 16 °nm^^, then gradually decreases to zero 
at higher angles. Allowing relaxation in the twisted nan- 
otube structure has a strong influence on the value of the 
band-gap; the gap increases linearly (but with a different 
slope) with twisting angle up to 14 °nm~^, then reaches 
a stable value ( « 0.3 eV) at higher angles. The first, 
linear, region corresponds to the opening of a gap pro- 
duced by the strong effects of symmetry breaking on the 
frontier tt and tt* orbitals. This initial linear increase ob- 
tained in constrained structures is in general agreement 
with the predictions of low energy theory |^ which, how- 
ever, does not take relaxation into account. For higher 
twisting angles, the band-gap decrease observed for the 
constrained geometry is mainly due to the presence of 
highly destabilized a — cr*-orbitals near the Fermi level 
that emerge from the increasing {sp)a — isp)a overlap. 
The plateau observed in relaxed structures is the result 
of the collapse of the tubular structure into a flattened 
helix shape which decreases the large repulsive overlap of 
strained NTs (Fig.2c,d). 

The opening of a band-gap in originally metallic nan- 



otubes leads to a drastic modification of the electrical 
properties of nanotubes upon twisting by even low an- 
gles. In order to investigate the conduction of electrons 
in twisted (and relaxed) (6, 6) nanotubes, we have con- 
nected both tube ends (with dangling bonds on the end 
carbons) to gold electrodes. Each electrode is composed 
of a layer of 22 gold atoms in a (111) crystalline ar- 
rangement. The distance between the NT end and the 
gold layer is 1.0 A. Such a bonding configuration min- 
imizes the contact-resistance and emphasizes the trans- 
port properties of the nanotube itself. The Hamiltonian 
(iJjvr) and overlap matrices (Snt) of equation (1) are 
determined using EHM for the system Gold-NT-Gold. 
The transmission function T{E), that represents the sum 
of the transmission probabilities over the contributing 
nanotube conduction channels, is then computed from 
equation (3). The T{E) spectra are expressed with re- 
spect to the Fermi level {Ep is defined as eV) of indi- 
vidual twisted nanotubes. 
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FIG. 3. Band-gap opening in a (6,6) nanotube as a function 
of the twist angle for a constrained geometry (open circles), 
and a relaxed geometry (closed circles). 

Figure 4 shows the transmission T{E) spectra of nan- 
otubes twisted by 0, 1.6, 4.1 and 8.1°nm~^ angles (main 
panel). All these deformed tubes lie within the linear 
regime of the band-gap variation with twist angle. The 
variation of T{E) at Ep with twisting angle is shown at 
the upper-left panel, while the corresponding change in 
resistance at zero bias (= 12.9 'kD,/T{Ep)) is given at the 
upper-right panel. The spectrum of the perfect structure 
(0°nm"^) shows a T{E) at the Fermi level of about 1.2, 
leading to a resistance higher than expected for ballistic 
transport (where T{E) = 2.0). This reduction in trans- 
mission is due to the presence of a finite contact resis- 
tance. The increasing T{E) on the high binding energy 
side is due to the opening of higher conduction channels. 
The asymmetry of the transmission function T(E) above 
and below Ep depends on the NT-metal electrode cou- 



Electrical and Mechanical Properties of Twisted Carbon Nanotubes 



4 



pling (Au-NT distance). A larger Au-NT distance leads 
to a higher T(E) above Ep , and vise versa. Since the 
NT-metal electrode distance is kept fixed in all calcula- 
tions, this effect of contact resistance does not influence 
the effects induced by tube twisting. 
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FIG. 4. Computed transmission function of a twisted (6,6) 
NT. The upper-left inset shows the logarithmic decay of 
T{Ef) with twisting angle, and the upper-right inset gives the 
variation of the resistance at zero bias for the same twisted 
NT. 

Twisting induces important changes in the T{Ef), a 
property that determines the linear conductance of the 
system. Already a twist of 4° nm^^ (or 0.5°scction^^), 
reduces the transmission by a factor of 14. It is interest- 
ing, however, that essentially no change is seen in T(Ep) 
at low twisting angles (< 2°nm^^) despite the fact that 
a finite band-gap (« 0.04 cV) is calculated for the free 
NT. We attribute this behavior to the metal-induced gap 
states (MIGS) produced by the interaction of the short 
(96 A) NT segment with the gold metal electrodes Q 
and the presence of charge-transfer doping. The MIGS 
tend to bridge the distortion-induced gap when the lat- 
ter is small [g^l- For twisting angles larger than 2°nm~^, 
the transmission decreases exponentially (see upper-left 
in Fig. 4). At 8°nm~^, T{Ep) has been decreased by 
more than 2 orders of magnitude, and if we consider 
transport into a strongly twisted, helical shape nanotube 
(see Fig. 2 [c,d]), the decrease in T{E) is about 5 orders of 
magnitude. A twisting-induced band-gap is, most likely, 
responsible for the recently documented field-effect tran- 
sistor action of a twisted, normally metallic multi-wall 
nanotube 

In contrast to our findings on axially bent nanotubes 
|0| where for a (6, 6) tube the decrease in transmis- 



sion occured at energies below Ep, and was traced to 
curvature- induced localized a-ir orbital mixing, the sharp 
decrease in T{E) upon twisting straddles Ep and is in- 
duced by symmetry breaking all along the NT. 

We thank Dr. Thomas Heine for providing us the tight- 
binding DFT code, and Dr. Frederic Lesage for very 
helpful discussions. 
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